Long-range entanglement generation via frequent measurements 
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A method is introduced whereby two non-interacting quantum subsystems, that each interact 
with a third subsystem, are entangled via repeated projective measurements of the state of the third 
subsystem. A variety of physical examples are presented. The method can be used to establish long 
range entanglement between distant parties in one parallel measurement step, thus obviating the 
need for entanglement swapping. 
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I. INTRODUCTION 

Entanglement is at the heart of quantum information 
processing 0] , and is a resource that allows quantum pro- 
cesses to outperform their classical counterparts for tasks 
such as^omputation Q , communication Q , and cryptog- 
raphy 4]. The standard way to generate entanglement 
between distinguishable and initially separable particles 
is to let them interact directly for a certain amount of 
time: any non-trivial two-body Hamiltonian is capable of 
generating entanglement in this manner . However, the 
direct interaction method poses limitations in the con- 
text of the generation of long-range entanglement, since 
in many systems the interaction strength typically de- 
creases at least as fast as some power of the distance be- 
tween particles. Here we introduce a different paradigm 
for entanglement generation: we show that it is possible 
to entangle two particles that never interact directly by 
means of repeated measurements of a third subsystem that 
interacts with both. In addition to its conceptual interest, 
we show that this scheme offers practical advantages for 
long-range entanglement generation. 

We remark that other alternatives to the direct in- 
teraction method for entanglement generation are well 
known. E.g., in Q linear optics is augmented by mea- 
surements to generate entanglement and perform quan- 
tum computation. Measurements of the phase of light 
transmitted through a cavity have also been proposed 
in Q as a method to prepare entangled states and im- 
plement quantum computation in the case of atoms in 
optical cavities. Decoherence can also be used to create 
entanglement: e.g., atoms in a leaky cavity can become 
entangled conditioned on null detection at a photo de- 
tector placed outside the cavity H, 0| ■ In Ref . |0 en- 
tanglement is generated in a similar manner conditioned 
instead upon detection. 

Here we present a rather general theoretical framework 
for measurement-generated entanglement. Our scheme 
is inspired by the recent work by Nakazato, Takazawa, 
and Yuasa (NTY) 1 1] , who investigated the effects of re- 
peated rapid projective measurements on one subsystem 
of a bipartite system. Provided that the resulting oper- 
ator on the unobserved subsystem has a non-degenerate 
largest eigenvalue, NTY showed that the unobserved sub- 



system is gradually projected into a pure state, indepen- 
dent of the initial state. This result leads to the question 
of whether a similar measurement scheme is capable of 
generating entanglement when the unobserved subsystem 
is itself multi-partite, non-interacting, and in an arbitrary 
initial state. We answer this question in the affirmative 
for a wide range of model systems, and establish general 
conditions for validity of the method. We note that an 
additional advantage of this measurement-based method 
for generating entanglement over the direct-interaction 
based method is that does not depend on sensitive tim- 
ing of interactions: instead, entanglement is gradually 
purified as the number of measurements increases, and 
(under appropriate conditions) can be made arbitrarily 
high. An important conclusion from our study is that 
the method can be used to establish long range entan- 
glement between distant particles, by measuring along a 
chain of intermediate particles. This provides an alterna- 
tive to entanglement swapping that does not scale with 
the chain length. A similar result - long range entangle- 
ment from local (bu t non-repeated) measurements - was 
obtained in Ref. [I^] for the ground state of an antiferro- 
magnetic spin chain. We now turn a detailed description 
of our method. 



II. PRELIMINARIES 

NTY considered the following scenario: Consider a bi- 
partite system composed of subsystems A, B, initially in 
the separable state pab{0) = \<l^){4'\(^ Pb, where A is in a 
pure state \(/)) and is arbitrary. Subsystem A is subject 
to projective measurements Pa — \4'){4>\ applied with 
period r. In between measurements the system evolves 
under the Hamiltonian 



H = Ha + Hb+ Hai 



(1) 



(resp., the sum of free Hamiltonians and an interaction). 
It can be shown that after M such measurements the 
state of subsystem B, given that all outcomes were 



M 



(2) 
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where Pm is the survival probabiUty, i.e., the probabihty 
of finding subsystem A in its initial state, 



PM = TrB[VB{T)''pBVUr)''] 



(3) 



where U — exp(— ir_ff) (we use units where h — 1), and 
where 



Vb{t) = (0|C/|0) 



(4) 



is an operator on subsystem B that is in general not 
Hermitian. However, one may still find left- and right- 
eigenvectors with complex eigenvalues, whose modulus 
can be shown to be bounded between and 1. The cen- 
tral result of NTY is that in the limit of large M and 
small but finite r, pb{Mt) tends to a pure state inde- 
pendent of B 's initial state. This result assumes that 
the largest eigenvalue A,„ax of Vb{t) is non-degenerate. 
The final pure state pb{Mt) is then the corresponding 
right eigenvector iMmax) and this outcome is found with 
probability 



Pi 



M 



\PB\Vn 



where li^max) is the corresponding left eigenvector. Note 
that because of r's finiteness the dynamics here is distinct 
from the Zeno effect and "quantum Zeno dynamics" |l3j| . 

We now extend NTY's model by allowing B itself 
to be composed of multiple non-interacting subsystems, 
B = {-Bi, -B2, -Bjv}, and pose the question: is the NTY 
measurement procedure capable of generating entangle- 
ment amongst _B's subsystems, assuming that they all 
interact with the "station" Al The situation is illus- 
trated in Fig. ^ To answer this questions we consider 
the following model: We set all internal Hamiltonians 
Ha = HBi — 0, so that H = Hab (equivalently, we 
can always transform to an interaction picture rotat- 
ing with the internal Hamiltonians; this will make Hab 
time dependent, but it will not introduce couplings be- 
tween subsystems Bi). We assume that the B subsys- 
tems are all qubits and the interaction between A and 
all the Bi is identical. We can then define total quasi- 
spin operators cr^ = X^il^i "'/j ^^'^ write the interaction as 
Hab = J2ae{x.v,z} ^aA^ ® a a {ha real) or, more simply, 



Hab = A- 



(6) 



where the parameters ha are included in the vector A . 
We do not restrict subsystem A. 



III. GENERAL THEORY 
A. A good basis 

We construct a basis that block-diagonalizes the effec- 
tive post-measurement evolution operator Vb{t). Let a 
basis for the A subsystem be the measurement state \4>) 




Measurement 



FIG. 1; Schematic illustration of our model: Systems Bi and 
B2 are coupled to system A; system A can be measured with 
a projective measurement. Note that no direct coupling exists 
between systems B\ and B2. 



and any set of states orthonormal to it, denoted 
I = l,...,d— 1. Let N be even and let a basis for the 
B subsystem be constructed out of the usual spin basis 
\S,Ms), where S is the total quasi-spin of the N parti- 
cles and Ms its projection along the z-axis. We denote 



(5) the (orthonormal) singlet states by \sj), j = 1, ...,D 



N 



7V!/[(iV/2-|- l)!(7V/2)!], and the remaining (orthonormal) 
states with 5 > by \tk), k = 1, ...,Kn, Kn = 2^-Dn. 

Proposition 1 : Consider the ordered basis 
{|(^),|</.^),...,|0j-_i)} ® {|si),...,|sD„),|ii),...,|tif„)}. 
In this basis we have the block-diagonal representation 



Vb{t) = 














(7) 



where Id^ is a 13 at x Ujv-dimcnsional identity matrix. 
The maximal eigenvalue of Vb{t) is 1, and is at least 
ZJjv-fold degenerate. 

Proof: The singlets states, \S,Ms) = |0,0), are anni- 
hilated by the ap when TV is even (since they are states 



with zero total quasi-spin). Therefore HAB\sj 
dependent of the state of subsystem A. 
basis Hab is thus represented as 



Hab = 



0, in- 
In the ordered 




(8) 



where the dimension of the upper-left block is x Da? 
and that of the lower-right block is Km x Km. Then 



U 



AB 



itHab _ 











U'ab 



(9) 



where U'j^^ = W'e~^^^W'^ is unitary and W is the uni- 
tary matrix that diagonalizes H'y^g-. W'^ H'j^^W' = A = 
diag(Ai, Aif„). Taking the expectation value wrt |0) 
we then find: 



D 

E 



VB{T) = y^\s,){s,\+Vi 



Bi 



(10) 



which is the claimed result, 
Y.k,k'=i\*k) {U'AB)^k.4,k' i^k'l a Km 



with V§ = 
X i^Tv-dimensional 



3 



matrix. That the maximal eigenvalue is 1 follows from 
unitarity of U'^^ jflij, and that it is at least I?jv-fold 
degenerate is immediate from the representation 0. 
QED 

Proposition 2 ■ Non-degeneracy of Vb (t) is a nec- 
essary condition for obtaining a pure state in the limit of 
large M. 

Proof: In the degenerate case it follows from Eq. Q 
that the method yields an equally weighted sum of degen- 
erate pure states corresponding to the maximum eigen- 
value of FbCt). QED 

Note that D2 = I. Hence in the totally symmetric case 
we have been considering so far, for N — 2 the method will 
generate a pure (maximally entangled) singlet state, pro- 
vided Vg has maximal eigenvalue with modulus smaller 
than 1. We give examples of corresponding Hamiltoni- 
ans below. On the other hand, D4 — 2 and hence due 
to degeneracy the method will not produce a pure en- 
tangled state for N > A. However, we can still generate 
pure-state entanglement in the N = 4 case by breaking 
the total permutation symmetry and only preserving the 
symmetry between 1,2 and 3,4 [e.g., having a coupling 
to subsystem A such that a/s = oi (cr^ -I- cr^) -I- 02 (crj^-fcr^), 
oi 7^ 02 7^ 0]. We can then project out a one-dimensional 
subspace |s)i2|s)34, \s)ij = (|0ilj)-|li0j))/V2 is the sin- 
glet state, which is entangled for 1,2 and 3, 4 but separa- 
ble across the 12 : 34 partition. Similarly, for = 2n > 6 
and ap — X^ILi '^ii'^p^"^ '^i 7^ "^i 7^ the projected 

state is |s)i2|s)34 |s)„-3,n-2|s)„-l,„. 

B. Invar iance 

The results of the above discussion are invariant as long 
as the Hamiltonians belong to the family {UsHABUy}, 
where Ub is an arbitrary unitary transformation of the 
Hamiltonian of subsystem B. This enables the method to 
generate entangled states equivalent under local transfor- 
mations. The invariance of HAB\sj) = is -ff^^jsj) = 0, 

where H'j^g = UbHabUb ^^"^ l^j) ^ t^ski)- Thus, e.g., 
in the case of TV = 2, we can generate the other three 
Bell states by applying Ub = ^2, ^2, ^2 (where X2 is the 
Pauli matrix Gx acting on subsystem i?2, etc.), which re- 
sults, respectively, in |s') ^ ;7|(|11) + |00)),^(|00) - 

|11)), + |10)) as the outcome of the method. 



C. Degenerate maximal eigenvalue 

What happens when there are other states, either sin- 
glets or coming from V^, that also have an eigenvalue 
with modulus 1? In this case NTY's result does not ap- 
ply and entanglement may or may not be spoiled (though 
it is always reduced from maximally entangled). Two 
examples will illustrate: (i) N — 2: Suppose there is a 



triplet state |00) that also has |Ai„ax| = 1- Then the 
resulting state is the mixture \s){s\ + 1 00) (001. This 
state is entangled as its partial transpose has a 
negative eigenvalue of —0.207. We will encounter this 
case in the Heisenberg model below, (ii) iV = 4: Sup- 
pose the two singlet states |si) — |s)i2|s)34 and |s2) = 
^[|i+)i2|t-)34 + |t-)i2|i+)34-2|to)i2|to)34] (where |t„) 
are triplets with projection quantum number a) appear 
with eigenvalue one, but no other states do. With respect 
to the 12 : 34 cut |si) is a product state but \s2) is clearly 
entangled. This state, just like in the previous example, 
has negative partial transpose, and we have entanglement 
across the 12 : 34 cut. These examples illustrate that de- 
generacy still allows for mixed state entanglement to be 
generated by our method. This is useful entanglement 
in the sense that it can be used for teleportation and all 
other QIP primitives [T5|. 



IV. EXAMPLES 

We now discuss examples, first limiting ourselves to A 
being a qubit and N — 2. Our task then reduces to calcu- 
lating the eigenvalues of the 3 x 3 matrix = {4>\U'ab\(I^) ■ 
However, this requires diagonalization of U'^g, an 6 x 6 
matrix, so cannot be done analytically in complete gen- 
erality. Our basis for A is {\(j>), |0^)} and for B is 

P^{\s) = i=(|01)-|10)),|L) = |00), 

\to) = -^(|01> + |lG»,|t+> = |ll)}}, (11) 

without loss of generality (recall the invariance discussion 
above). 

A. Axial symmetry model 

Suppose Az = and Ax,Ay satisfy [^3;,^^] = 
= (or Ax or Ay — 0), whence the Hamilto- 
nian reads H = [Ax{Xi + X2) + Ay{Yi + Y2)]/2. We 
first consider as a special case the XY model; Ax = JX 
and Ay = JY, relevant for quantum information process- 
ing in solid state 0, O and atomic [l^ systems. 
It follows after some algebra that = \A\^ H, where 
1^41^ = Al + Ay. Therefore, the evolution is 

Uab = e-^ = / - 2 ' sin^ ^ " ^^^^ l^" ' 

(12) 

which means that Vg is a function only of {(f>\H\(l)) and 
{(f)\H'^\(t)) . In particular, for the XY model we find, using 
\A\ = V2J, 



4 



Vb{t) 



tJ X1X2 + YxYi - {Z){Zi +Z2) .2 rJ 



V2 



V2 



I — sm \JLtJ 

2V2 



r 



X2) + JZ(Z\ + Z2) (A is a tunneling parameter and J is 
the Josephson coupling). Thus A = 2(A, 0, JZ), and we 
may use the general result (I12|l . replacing Ay by A^. If 
we now choose to be a Uz-eigenstate then Vb turns out 
to be unitary and hence all its eigenvalues have modulus 
one, and no entanglement is generated. However, if we 
choose 10) to be a (Ta;-eigenstate, we find 



Vb 



cos 
i 



sin^ 0(cos2 9X1X2 + sin^ 9Z1Z2) 



sin2(?!)cos6'(Xi + X2) 



FIG. 2; Concurrence (A) and success probability (■) Pm for 
l<^> = 75(10) + |1)) and rJ = 7r/2. 



where (X) = {(j)\X\(j)), etc. It is simple to check that, 
as required from our general result, Vb\s) = |s), i.e., the 
singlet state has eigenvalue 1. Whether this eigenvalue is 
degenerate is now seen to depend on the measurement of 
the A subsystem: If we choose to be a (Tz-eigenstate 
then the additional eigenvalues of Vb are found to be 
{1, cos \/2t J, cos V2tJ}, while if we choose |0) to be a 
or (Ty-eigenstate then we find that the additional eigen- 
values are {cos^ g, 1 + 3cos2g ± \/ cos 4g + sin'' q — 1}, 
where 9 = Thus, if we measure along z then there 
is never pure state entanglement, but at times r other 
than mr/V^J we have a mixed entangled state, since the 
eigenstate corresponding to the additional eigenvalue 1 is 
found to be |11) (and recall the discussion above). In the 
case of measurement along x or y the eigenvalues period- 
ically have modulus 1, and with the exception of those 
times degenerary is avoided and we do obtain a pure en- 
tangled state (in particular, it is simple to show that this 
is true for all r < \/2arccos(l/3)/ J). 

Next we consider the performance of the scheme after a 
finite number of steps (Fig. To do so we calculate the 
concurrence [23| after M steps of evolution/measurement 
for Jt = 7r/2 < -\/2arccos(l/3), in the case of |0) a Ux 
eigenstate. In this case, as shown above, the singlet state 
is generated by the scheme. We also plot the survival 
probability Pm as given by Eq. Q . After as few as three 
measurements the concurrence is essentially unity, indi- 
cating that the system is already in the desired singlet 
state. The scheme does, however, not work with unit 
probability: The probability Pm of having projected sys- 
tem A on the desired ax eigenstate converges to ^ 0.2 
for the values above, a value reached for M — i. 

As another special case we consider a simplified form 
of the d-wave grain boundary qubit ^21|: H = A(Xi -|- 



where (j) = rV A^ + ,P and tan 6* = J/ A. It is again 
easy to check that Vb\s) = \s). One of the other 
three eigenvalues is 1 — 2 sin^ cj) sin^ 9 . The other two 
are complex conjugates and for short times t such that 
tan^ (/) < 4 cos^ 6*/ sin'' have the same amplitude (1 — 
2sin^ 0sin^ 6*)^^/^. For instance, when 9 = 7r/4, the am- 
plitudes of the three eigenvalues are cos^ cj), cos 4>, cos (j>, 
which vanish rapidly when raised to the power of the 
number of measurements. 

Finally, note that the results presented so far are not 
restricted to A being a qubit: the unitary operator H12|) 
is valid even when ^ is a multi-level system, as long as 
the condition [A^,,^^] = [Ay,Al] = is satisfied. 



B. Heisenberg model 

We now consider the Heisenberg interaction H ~ 
J~a ■ ("(T 1 + ~G2)- After showing that {H + J~a \ ■ "ti 2)^ — 
9 and \H, J'a i ■ 'a 2] = 0, it is simple to compute 
^-it{h+j u 1- G 2) ^ from which we directly obtain 



Uab = e-""" =e'^'"'"^{cos{-iTj) 

-i sin(3Tj)["(T • ("(T 1 + ~a2) + ■ 0^21/3}. 

Vb then is identical to Uab provided one replaces with 
(01 "ct 10) the operator 'a in Uab- The four eigenval- 
ues are {1, e"^''^'', e''^'' (cos 3r J ± zi sin3T j)}, the last 
two having the same amplitude 1 — | sin^ 3tJ < 1. The 

corresponding eigenvectors are, respectively, (3 [Eq. 
Since the first two eigenvalues both have magnitude 1, we 
have the case discussed above: a pure entangled state can- 
not be projected out by our method, but we can prepare a 
mixed entangled state that is useful for all QIP protocols. 
Note further that these Heisenberg model results hold for 
the entire class of Hamiltonians H — A^, ■ ilri + ^2), 
when is generated from an arbitrary 2-dimensional 
unitary transformation LOi, A^, = UaJ '^U\ (this is 



5 



different from the invariance under a rotation of the B 
system considered above). 

C. Bosonic media 

We now consider a photon or phonon (A) interacting 
symmetrically with two identical qubits (B). The Jaynes- 
Cummings Hamiltonian is 

H = eb^b + g{al + a|) + i[b{a+ + a+) + b\a^ + a^")] 

where 6 is a bosonic annihilation operator, and whence 
A = {J{b + b^),iJ{b — b^),g)- Such a model can easily 
be realized using microwave cavity QED 22] (two atoms 
in one or two cavities) or trapped ions [23 • We can 
exactly diagonalize H by writing it as a direct sum over 
three-dimensional matrices in the basis |n — 1, t+) , |n, to) 
and \n + l,t_), where \n) are number states and \ta) are 
triplet states of the two qubits. Note that J\f = [b^b + 
{af + ct|)] is a conserved quantity and hence H is block 
diagonal in its eigenvalues. For simplicity we consider 
the case with g = e and measure the single photon state 
{(j)) — b^ |0) — |1) (this projects onto a single block of 
the infinite matrix Uab)- In this case we readily find the 
already diagonal 

Vb{t) = diag{l,e ^""^ 

5 

e-*" cos(%/6Tj),cos(%/2rJ)} 

in the ordered basis {|l)|s), |1, <+) , |1, Iq) , |1, L)}. Thus, 
as long as we make sure that r < 2tt/VT0J the method 
will project out a pure singlet state. It follows from our 
discussion of invariance above that another Bell state 
-^(|01) + |10)) can be projected out if the two qubits 
couple to the photon or phonon with opposite signs. 

D. Multilevel Systems 

Let us now consider a rather general, though abstract 
multilevel case. Assume that subsystem B consists of 
two particles that have M > 2 levels each, and that the 
interaction Hamiltonian with subsystem A (also an M- 
level system) is of the form 

M 

H=Y,A,,{0l^+0ll 

where Oij — Eij — Eji and Eij = \i) {j\ is a matrix whose 
elements are zero everywhere except for a 1 at position 
(i,j). Namely, there is an SO{M) symmetry for odd M 
or an Sp{M) symmetry for even M . In this case H has a 
single non-degenerate eigenvector with zero eigenvalue. 
This state is entangled. For instance, when = 3, we 



have two qutrits, which can be represented in the spher- 
ical basis {| — 1) , |0) , The state with zero eigen- 
value is 

|vl') = (|l)j-l),-|0)j0)2 + 1-1)1 |1)2)/V3, 

and is maximally entangled. In general, the one- 
dimensional subspace containing the state with zero 
eigenvalue is the irreducible representation (0, 0, 0) of 
SO{M) or Sp{M). This state is the generalization of the 
singlet state that arose in the general theory and exam- 
ples treated above when the subsystems were qubits. 



V. PERTURBATIVE TREATMENT 

It is apparent from the above examples that the suc- 
cess of our method depends on keeping the period t 
between two measurements short. For instance, in the 
bosonic example, we require r < 2t:/^/\QJ in order to 
prevent the appearance of another eigenvalue with am- 
plitude one. Let us therefore consider a short-time ex- 
pansion of Vb (t) . To first order 

Vb{t)^' - {I-lT{cf>\HAB\<P)f' 

exp(-ii (01 Hab 10)), 

where Mr = t (constant), i.e., the evolution is unitary. 
This is the quantum Zeno effect, which completely de- 
couples the interaction among all parties of subsystem B, 
so that no entanglement can be generated in this limit. 
Thus the dominant contribution to entanglement gener- 
ation originates from the second-order term, which con- 
tains the self- correlation of subsystem B. Consider for 
simplicity the case (01 Hab |0) = (as in our last exam- 
ple). Letting Mt'^/2 = (constant), we have 

Vb{tY'^ (,-'t^{'t>\HlB\<t>) 

All eigenstates of (0| |0) with nonzero eigenvalues 
are rapidly suppresed as t (in practice M) increases, 
while those with eigenvalue zero survive. Since it is much 
simpler to analytically calculate (0| |0) than Vb{t), 
this perturbative method provides a relatively simple 
tool for estimating the possibility of entanglement gen- 
eration via our method, for complicated systems. We 
consider spin-orbital coupling as another illustrative ex- 
ample: H — hLx{Xi + X2) + hLy{Yi -\- 12), denoting two 
spin qubits that couple with the same orbital angular 
momentum. If we measure the eigenstate |0) = |/, 0) 
of the orbital component, we find (0| |0) = and 
(01 i72 ^ 1^1(^1 + i)(2 + X1A2 + FiFa)- The eigenval- 
ues are hl{l + l){0, 2, 4, 2} with respective eigenvectors /?, 
meaning that the singlet state \s) is projected out in the 
second-order analysis. Note that this example is qualita- 
tively different from the previous one since subsystem A 
here does not refer to a physically separate particle. 
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VI. PREPARATION OF LONG-DISTANCE 
ENTANGLEMENT 

Finally we come to our main result: the generation 
of long-distance bi-partite entanglement. Since in our 
discussion above there was no restriction on the size of 
the A subsystem, in a sense long-distance entanglement 



generation already follows from the results above. How- 
ever, it is important to specify how the measurements on 
A can be carried out. Assume that the A subsystem is 
composed of — 2 qubits 2, 3, — 1 and we wish to 
entangle the two _B-subsystem qubits 1, A^. We consider 
again as an illustrative example a simplified form of the 
(i-wave grain boundary qubit |2ll | : 



Uab = exp[-i Jr(Xi + Xn) - iAr(ZiZ2 + Z2Z3 + hs + Zn-2Zn-i + Zm-iZm)] 



where /is — X^i^s^ -^j-^j+i- We choose \4>) to be the state |+)2 \R) where \R) = \+)^ l+)4--- \+)n-2^ where 

1+) = (|0) + |l))/%/2. It follows after some calculations that 

Vb = i{a[cos^ (/) - sin^ (/)(cos^ ^Xi^AT + sin^ 6'ZiZjv) - «^sin20cos0(Xi + Xjv) 
-|-5[cos^ (f) - sin^ (/)(cos^ OXiX^ - sin^ BZiZn) - sin20cos6'(Xi + Xn)} 



where 



(and are real numbers for odd N), (f> — rvJ^ + A^, 
tan 6' = A/ J. Two of the eigenvalues and eigenstates 
are: 

El = i{a + 6(l-2sin2</>sin2 0)};-^(|liOjv> - lOiljv)), 

E2 = i{& + a(l-2sin2 0sin2 0)};-^(|OiOjv> - lliliv)). 

Thus, as desired, entanglement is generated between the 
B subsystem particles 1, N. Note that since the measure- 
ment of all the particles of the A subsystem is carried out 
in parallel this method for long-range entanglement gen- 
eration is independent of the distance N. 

The amplitudes of the other two eigenvalues are 
\/EiE2 in the case of odd N > 5. Since max(|i?2|, > 
y/EiE2 for sufficiently short time, and Ei — E2 = 
2 (sin (/))^ (sin 6*)^ (a — 5) as long as a ^ 5, the proce- 
dure prepares either --i=(|liOAr) — IOiIat)) or --i=(|0i07v) — 
IIiIat)). We find, e.g, for TV = 5: a = cos2At, b = 1 
and --i=(|00) — |11)) is projected out, while for N = 7: 

a = i(3 -I- cos4Ar), b = cos2At and ^(|01) - |10)) is 
projected out. The case of even N is more complicated 
since the amplitudes of other two eigenvalues depend on 
the values of a and b. E.g., in the N = 4, Vb has the 



same form as above but a — e^'^"^ and b = e*'^'^; all 
roots have the same amplitude and no pure state will be 
projected out. 

VII. CONCLUSIONS 

We have introduced a measurement-based method for 
entangling two systems that only interact indirectly, via 
a third system. Applications range from solid state to 
atomic quantum information processing. One of the ad- 
vantages of this measurement-based method of entangle- 
ment generation is that it does not depend on precisely 
timed interactions, as are interactions-based schemes. 
The method can be used to prepare arbitrarily long- 
distance entanglement via a chain of intermediate par- 
ticle in one parallel measurement step. This may have 
useful applications in reducing the latency overhead in 
quantum communication in quantum computing archi- 
tectures. An interesting open question which we leave for 
future research is whether the same method can be used 
to apply quantum logic gates between non-interacting 
particles. 
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